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摘 要: 使 用 开源 流体 动力 学 程序 库 OpenFOAM 对 竖 直 圆 管内 超 临界 压力 RP-3 对 流传 热 进行 数值 模拟 ， 研 究 重 点 为 起 始 加 
热 段 的 传 热 恶 化 现象 。 本 文选 择 k-e-v-f F k-w SST 两 种 低 Re 数 汕 流 模型 进行 计算 并 与 实验 数据 进行 对 比 ， 讨 论 了 两 种 模型 
对 于 高 壁面 热流 条 件 下 传 热 恶 化 现象 的 预测 能 力 。 通 过 对 流 场 及 温度 场 的 分 析 对 起 始 加 热 段 传 热 恶 化 现象 给 出 机 理解 释 ， 最 
后 考察 了 变 物性 及 浮 升力 的 影响 。 结 果 表 明 : k-o SST 模型 无 法 模拟 传 热 恶化 ，k-e-w-f 模 型 能 够 定性 描述 起 始 加 热 段 内 壁 温 
的 迅速 升 高 。 由 流体 密度 变化 及 重力 引起 的 浮 升力 作用 是 导致 起 始 加 热 段 传 热 恶化 的 重要 原因 。 
关键 词 : 超 临界 压力 流体 ， 对 流 换 热 ， 数值 计算 ，RP-3 
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Abstract: The open source CFD toolbox OpenFOAM was used to simulate convective heat transfer of RP-3 flowing 
in a vertical tube under supercritical pressure. Heat transfer deterioration (HTD) in the initial heating region was 
studied. Two types of low-Reynolds number turbulence models (k-z-v^-£ and k-co SST model) were used in the 
simulation. Numerical results were compared against the experimental data to examine the accuracy of turbulence 
models in predicting HTD under large wall heat fluxes. The flow field and temperature field were analyzed to 
propose a mechanism for explaining HTD. The effects of property variations and buoyancy forces were also 
investigated. It was concluded that the k&-co SST model was incapable to predict HTD. The &-e-v^-f model was able to 
qualitatively describe the rapid increase of inner wall temperature in the initial heating region. The buoyancy forces 
caused by density variation and gravity is the main cause of HTD. 
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超 燃 冲 压 发 动机 是 吸 气 式 高 超声 速 飞行 器 的 室外 壁 训 米 量 级 的 冷却 通道 吸收 壁面 热量 ， 再 喷 
核心 动力 部 件 ， 长 时 间 、 高 强度 的 超声 速 燃 烧 与 入 燃烧 室 进 行 燃烧 。 为 了 提升 换 热效率 并 避免 相 
严重 的 气动 加 热 会 带 来 极 高 的 热 负 荷 目 20. ETE 变 的 发 生 ， 冷 却 通道 内 的 压力 高 于 燃料 的 临界 压 
赫 数 飞行 速度 下 ， 来 流 空气 的 滞 止 温度 也 随 之 升 力 加 。 超 临界 压力 流体 传 热 的 特殊 性 在 于 其 热 物 
高 ， 因 而 采用 空气 作为 冷却 介质 的 冷却 方式 无 法 性 会 随 压力 和 温度 的 变化 而 发 生 很 大 的 改变 。 相 
应 用 于 超 燃 冲压 发 动机 的 热 防护 中 四 。 目 前 ， 采 关 研究 通常 将 一 定 压力 下 定 压 比 热 〈Cp) 峰值 所 
用 机 载 碳 氧 燃料 作为 冷却 剂 的 再 生冷 却 技 术 被 认 对 应 的 温度 定义 为 拟 临界 温度 〈7he) o EAR 
为 是 实现 超 燃 冲压 发 动机 主动 热 防护 最 为 有 效 的 界 温度 附近 很 小 的 温度 范围 内 ， 流 体 的 热力 学 及 
技术 手段 之 一 外。 输 运 性 质 《〈 密 度 、 粘 度 、 定 压 比 热 、 导 热 系 数 


在 再 生冷 却 过程 中 ， 碳 氧 燃料 首先 流 经 燃烧 等 ) 会 发 生 剧烈 的 变化 。 此 外 ， 由 于 流体 密度 变 


化 导致 的 浮 升力 及 热 加 速效 应 也 会 对 流 场 造成 影 300 _ mm， 加热 段 上 、 下 游 分 别 设置 90 mm 和 60 
响 。 以 上 因素 综合 作用 ， 使 得 超 临 界 压 力 流体 的 mm 的 绝热 段 。 实 验 采 用 对 不 锈 钢管 直接 通电 的 
对 流 换 热 过 程 十 分 复杂 [9。 方式 进行 加 热 以 实现 近似 均匀 壁面 热流 条 件 ， 更 
为 了 掌握 超 临界 压力 下 碳 氢 燃 料 的 对 流传 热 多 信息 详 见 文献 00。 本 文 数值 模 拟 工 况 参数 总 结 
特性 ， 国 内 外 学 者 已 开展 了 大 量 研究 工作 。Zhang 如 表 1 所 示 。 
AD. 1 对 超 临界 压力 RP-3 航空 煤油 在 光滑 圆 管内 1.2 RP-3 热 物性 
竖 直 下 降 流 的 传 热 特性 进行 了 实验 研究 ， 将 传 热 RP-3 航空 煤油 的 临界 点 为 03: T=645.04 K, 
过 程 分 为 四 个 阶段 : 起 始 加 热 区 、 正 常 换 热 区 、 p.-2.33 MPa. RP-3 航空 煤油 热 物性 与 输 运 性 质 的 
传 热 强化 区 和 传 热 恶化 区 ， 并 基于 实验 数据 拟 合 选取 依据 与 文献 0 相同 。 其 中 密度 、 动 力 粘度 和 
得 到 了 传 热 实验 关联 式 。Li 等 上 0 对 超 临 界 压力 定 压 比 热 采 用 Zhang 等 4648 采用 实验 手段 测定 的 
RP-3 航空 煤油 在 光滑 圆 管内 竖 直 上 升 流 的 对 流传 物性 数据 。 导 热 系 数 采用 Zhong 等 外 提出 的 10 组 
热 进行 实验 研究 ， 讨 论 了 质量 流 率 、 热 流 密度 、 分 替代 模型 依据 扩展 对 比 态 原理 计算 得 到 。 


压力 以 及 进口 温度 对 于 换 热 的 影响 规律 。 刘 波 等 
021 对 超 临界 压力 下 正 奖 烷 在 竖 直 圆 管内 的 对 流传 


表 1 数值 模拟 工 况 参数 


Table 1 Parameters of simulation 


热 进行 实验 研究 ， 讨 论 了 两 种 管 径 条 件 下 《内 径 


0.95mm 和 2mm) 浮 升 力 及 热 加 速 对 于 换 热 的 影 T G 

响 程度 与 作用 规律 。Zhong 58051 对 超 临界 压力 下 e NES i M E a a 

RP-3 航空 煤油 在 水 平 圆 管内 的 对 流传 热 进行 了 数 3 300 45 — 15727 

值 研究 ， 发 现在 较 高 壁面 热流 密度 条 件 下 ， 当 辟 309 35 — 1523 

温 略 高 于 拟 临界 温度 时 会 出 现 传 热 恶 化 现象 ， 其 MEE. EE 

JAN 3r RE TR Dot Pts fe 4) ARK. Jiang 等 pr 2 550 373 — 15727 
C2 04 对 超 临界 压力 下 JP-8 航空 煤油 在 竖 直 圆 管内 
一 的 对 流传 热 进行 了 数值 研究 ， 研 究 讨论 了 不 同 油 lop me | Loxo? 4020 
C2 流 模型 (标准 ke 模型 使 用 标准 壁面 函数 、RNG de —— Iri sei heat capacity | 9^ | 018 

k-e 模型 使 用 增强 型 壁面 函数 以 及 大 SST 模型 ) ee a 

对 于 传 热 恶化 现象 的 预测 能 力 。 研 究 认为 粘度 随 ut ms. 
^ 温度 的 剧烈 变化 是 高 热流 密度 工 况 下 传 热 恶化 发 - Uu" 
Zoo 生 的 首要 原因 。 Be om ¥ 
: 本 文 对 超 临界 压力 下 RP-3 航空 煤油 在 竖 直 光 "T pq 

滑 圆 管内 的 对 流传 热 进行 数值 研究 ， 重 点 关注 起 T p " "EL 

始 加 热 段 的 传 热 恶化 现象 。 文 中 对 两 种 低 Re 数 7". Soo 49 500 pe NES $00" od 

RANS MRE Ck-e-v2-f fU k-o SST) 对 于 起 始 加 

热 段 传 热 恶化 现象 的 预测 能 力 进行 研究 ;通过 对 图 1 RP AVE (g-5MPa) 

流 场 与 温度 场 的 分 析 对 传 热 恶人 化 给 出 机 理解 释 ; Fig. 1 Thermophysical properties of RP-3 (p-5MPa) 

最 后 讨论 了 变 物性 及 浮 升 力 对 于 起 始 加 热 段 传 热 2 数学 模型 

特性 的 影响 。 

2.1 控制 方程 
1 物理 模型 航空 煤油 管内 对 流传 热 是 典型 的 耦合 换 热 问 
ii EENGDEE 题 ， 本 文采 用 分 区 求解 、 边 界 看 合 的 方式 进行 处 
本 文 依据 Zhang 等 00 的 实验 结果 设置 数值 模 HH 由 0 

拟 工 况 ， 对 超 临界 压力 RP-3 航空 煤油 起 始 加 热 自 Bee | 

的 传 热 特性 进行 研究 。 实 验 中 ， 航 空 煤油 自 上 向 固体 域内 热传导 的 控制 方程 为 傅立叶 定律 : 

下 流 经 竖 直 放置 的 光滑 圆 管 (内径 1.805 mm， 外 actui (D 

1$ 22 mm)。 实 验 管 总 长 450 mm， 其 中 加 热 段 长 Ox, 


式 中 ,4 为 固体 导热 系数 ，$ 为 内 热源 。 
流体 域 控 


EG ))20 2 
a (ou) (2) 
动量 方程 
6 EN ôu, 
ar, P") ar ox Talt 
(3) 
能 量 方程 
H A ôH 
ap H - 
a 0n ) = 2 Nor 
(4) 


以 上 各 式 中 ，p 为 密度 ,4 与 jy 分别 为 分 子 粘 
度 与 湛 流 粘度 ，g 为 重力 加 速度 ; H KRAE, 
Pr 5 Pr 分 别 为 分 子 普 朗 特 数 和 满 流 普 朗 特 数 ， 

本 文中 Pr 二 0.85。 能 量 方程 中 忽略 了 重力 作 功 和 
粘性 耗 散 项 。 各 汕 流 模型 中 Am 的 定义 将 在 下 一 节 
中 给 出 。 

2.2 imi! 

超 临界 压力 流体 管内 对 流 换 热 与 近 壁 区 流 场 
密切 相关 ， 获 得 精确 的 近 壁 区 流 场 对 于 相关 研究 
非常 重要 。 因 此 ， 选 择 合理 的 湛 流 模型 对 于 超 临 
界 压力 流体 传 热 数值 研究 来 讲 是 十 分 重要 的 。 本 
文选 择 两 种 在 超 临 界 流体 传 热 领域 应 用 较为 广泛 
的 低 Re 数 RANS 江 流 模型 ，k-e-v?-f 模 型 I 和 k-w 
SST 模型 EI。k-e-y?-f 模 型 形式 如 下 : 

大 输 运 方程 : 
Opuk) ô M 


Ok 
十 +P +G, -pE 
d x; [Qu a ax k TOTP 


i 


(5) 

€ 输 运 方程 : 
O(puié) E CaP, *G,)- C. pe $ 6 E Ht 0€ 
Ox, T, Ox, c, Ox, 


i S i 


AP, PuZyiisz Be XJ: Gy 为 浮 升力 生成 项 ， 
采用 广义 梯度 扩散 假设 (Generalized gradient 
diffusion hypothesis, GGDH) 构建 C0， 其 形式 分 
别 如 下 : 


P=-puu, —- (7) 


ad pluh D (9) 
Oy, e, k 
n rp i 
poe £746 ye ‘i e — se (105 
EE A 


u, = pC, YT, (11) 


式 中 ， 和 常数 双 与 到 的 形式 详 见 参考 文献 09]。 


k-o SST 模型 形式 如 下 : 
大 输 运 方程 : 
O(pu,k) -p H 

Ox, 


co 输 运 方程 : 
O(pui®) 
Ox 


lure) ]-f'pko (12) 


0 Oo 
= aps? = Bpa’ t —[(u F O pkh ) —] 
， ôx, Ox, 
i-is A (I3 
"ra x, 


式 中 ， B nmin(p 105 pka) 


P,—u,— T). 
ey 


pl He (14) 
MEE max(d,o, SF,) 


k-o SST 模型 中 的 常数 作 按 下 式 计算 : 
$-Ó4F *9-RF) (15) 
混合 函数 F 定义 如 下 ; 


式 中 ，CD，= max| 2po,, 10"|; d 5 
Q9 Ox, Ox, 


"UE 500v. 4p0,,k 


Bey yo CD,y 


o2 


F = tanh |o 


C16) 


从 TIIK k-o 模型 和 标准 k-e 模型 的 常数 ， 详 见 


参考 文献 Bol。 
2.3 数值 方法 

数值 计算 使 用 开源 计算 流体 力学 程序 库 
OpenFOAM (版 本 号 2.4.0.) 对 偏 微分 方程 组 进行 
求解 。 控 制 方程 采用 有 限 体积 法 离散 。 动 量 方 
程 、 能 量 方程 中 对 流 项 采用 QUICK 格式 ， 其 余 各 
输 运 方程 中 的 对 流 项 采用 一 阶 迎 风格 式 。 采 用 
SIMPLE 算法 处 理 压 力 与 速度 耦合 。 当 各 求解 量 


残 差 小 于 105 时 认为 计算 收敛 。 

计算 域 尺 寸 与 物理 模型 相 一致 ， 加 热 段 上 、 
下 游 均 设置 绝热 段 以 保证 进入 加 热 段 时 流动 充分 
发 展 以 及 数值 解 的 稳定 。 为 保证 近 壁 区 流 场 的 精 


确 求解 ， 对 流体 域 壁面 附近 处 网 格 进行 加 密 。 本 s 
文 所 有 计算 工 况 均 满足 壁面 附近 第 一 个 节点 处 540 Y Zhang et al. (2012) 
mesh number 72000 
y*<0.5. 520 - - - mesh number 108000 
计算 边界 条 件 设 置 如 下 : 固体 域内 设置 均匀 一 一 mesh number 216000 
体积 热源 以 模拟 均匀 壁面 热流 条 件 ， 管 外 壁 与 两 3000. 20 40 60 80 100 120 140 160 180 
端 均 设置 为 绝热 。 流 体 域 入 口 处 设置 均匀 速度 、 x/d 
温度 及 各 消 流 输 运 量 ， 设 置 轴 向 压力 梯度 为 零 ; 图 2 不同 网 格 计算 内 辟 温 与 实验 数据 对 比 
出 口 处 设置 各 输 运 量 轴 向 梯度 为 零 ， 恒 定 压力 。 Fig.2 Comparison of calculated inner wall temperature 
壁面 处 各 输 运 量 边界 条 件 如 表 2 所 示 。 against the experimental data using different meshes 
> Mp WEIT 关 性 分 uk] 
m c T 3 计算 结果 与 分 析 
唱 型 对 工 号 工 况 进行 计算 ， 三 种 网 格 划分 方 
Vo 案 总 网 格 数 分 别 为 72000、108000 及 216000。 计 3.1 不 同 淇 流 模型 对 起 始 加 热 段 传 热 恶 化 的 预测 
) 算 工 况 参 数 如 下 : 压力 5 MPa， 质 量 流量 1572.7 图 3 给 出 ILIV 号 工 况 下 采用 kw SST 模型 和 
kg/m2.s， 入 口 温 度 473 K， 热 流 密度 300 kW/m?, ke-V?-f 模 型 计算 得 到 的 内 壁 温 与 实验 数据 的 对 
竖 直 下 降 流动 。 图 2 给 出 3 套 网 格 计算 结果 与 实 比 。 下 文 各 图 中 ，x/4d=0 表示 加 热 段 的 起 始点 。 模 
验 数 据 的 对 比 。 该 工 况 处 于 涨 流 区 主流 雷诺 数 拟 工 况 流 体 进 口 温度 较 低 ， 密 度 较 大 ， 流 速 较 
范围 14900-45000)， 内 壁 温 沿 流动 方向 单调 上 低 ， 管 入 口 流 动 属于 过 渡 区 进口 雷诺 数 约 为 
升 。 计 算 结果 表明 : 当 网 格 数 增加 至 108000 时 已 6600)。 实 验 数 据 显 示 ， 在 管 加 热 段 起 始点 附近 内 
经 可 以 满足 计算 结果 与 网 格 数量 无 关 ， 因 而 采用 壁 温 迅速 上 升 出 现 局 部 峰值 ， 发 生 传 热 恶化 。 在 
— ”该 方案 。 计算 内 壁 温 与 实验 数据 符合 良好 ， 也 验 换 热 强度 恢复 后 ， 内 壁 温 开 始 下 降 ， 再 沿 流动 方 
re 证 了 本 文 模型 与 数值 方法 的 正确 性 。 向 单调 上 升 。 如 图 3 所 示 ， 采 用 k-w SST 模型 计 
à —M—— — 算得 到 的 内 壁 温 治 流动 方向 单调 上 升 ， 未 观察 到 
局 部 峰值 ， 无 法 描述 传 热 恶 化 现象 。 采 用 k-e-v-f 


Table 2 Wall boundary conditions in the fluid 


模型 计算 得 到 的 内 壁 温 与 实验 数据 吻合 良好 : 内 


pem region 


"e 壁 温 沿 流动 方向 首先 迅速 上 升 ， 随 后 略微 下 降 ， 
a = k-e-v?-f 模型 k-o SST 模型 S ANT FE RO 
! NC ——"—Ó 再 沿 流动 方向 单调 上 升 。 可 见 ， 该 模型 能 够 定性 
ge ce id din 描述 起 始 加 热 段 内 壁 温 的 变化 趋势 。 
Pan eds k-o SST TUS Bri Vise ee np EERI, Tf k-e- 
zeroGradient zeroGradient m TREA 5 po HF Hy 
vf SCRI f FE] v" 作为 速度 尺度 以 正确 描述 壁面 处 


k fixedValue le-10 fixedValue le-10 


满 流 的 衰减 ， 部 分 考虑 了 潮流 各 向 异性 的 影响 。 


E zeroGradient / A : r HEA 

在 本 文 所 涉及 的 问题 中 ，k-e-w-f 模 型 获得 了 更 好 
[o / fixedValue le-10 $ 

的 结果 。 
y? fixedValue le-10 / 


f fixedValue le-10 / 
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Fig.3 Comparison of calculated inner wall temperature 


using different turbulence models against the experimental 
data 

图 4 进一步 给 出 采用 k-e-v f 模型 计算 得 到 的 
起 始 加 热 段 对 流 换 热 系数 与 实验 数据 的 对 比 。 由 
图 可 知 ，k-e-w-f 模 型 能 够 良好 地 模拟 起 始 加 热 段 
传 热 恶化 现象 ， 下 文 将 基于 该 模型 的 计算 结果 对 
传 热 恶 化 进行 分 析 。 
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Fig.4 Comparison of calculated heat transfer coefficients 


using the k-e-v?-f model against the experimental data 


3.2 起 始 加热 段 传 热 恶化 机 理 


图 5 给 出 本 号 工 况 下 采用 ke-v?-f 模 型 计算 得 


到 的 起 始 加 热 段 内 壁 温 及 对 流 换 热 系数 与 实验 数 
据 的 对 比 。 为 便于 分 析 , 沿 流动 方向 选取 4 个 截面 
MLE: P1-P4 CPlu/d-8,  P2:x/d-15,  P3:x/d-23, 
P4:x/d=40)。 其 中 ，P1-P3 位 于 传 热 恶 化 区 ，P4 位 
于 传 热 恶化 


HTC/W-m^- K^ 
N [99] 
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© © 
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后 换 热 的 恢复 区 。 
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Fig.5 Comparison of calculated inner wall temperature 


and heat transfer coefficient using the k-e-v?-f model against 


图 6 给 出 P1-P4 截面 处 流体 轴 
与 温度 〈 图 b) 的 分 布 情况 。 从 图 6 (a 可 以 看 
出 ， 在 传 热 恶化 区 (P1-P3)， 热 边界 层 内 流体 轴 
向 速度 呈 抛 4 


the experimental data (Case II) 


可 速度 (图 a) 


鸭 形 分 布 。 对 比 观察 图 6(b)， 此 时 壁 


面 与 主流 流体 温差 较 大 ， 主 流 区 流体 温度 变化 很 
小 。 其 原因 是 此 时 热 边 界 层 内 流动 呈 层 流 状态 ， 
热量 难以 从 壁面 传递 给 主流 区 流体 。 此 外 ， 热 边 
慨 的 厚度 治 流动 方向 逐渐 增加 ， 传 热 热 阻 不 断 
增 大 ， 从 而 导致 了 传 热 恶化 的 发 生 。 随 着 流体 温 


界 


界 


度 的 升 高 ， 密 度 减 小 ， 流 体 速 度 不 断 增 加 ， 热 边 
层 内 流动 开始 向 消 流 状态 过 渡 〈P4)。 由 图 6 


(b) 可 知 ，P4 处 壁面 与 主流 流体 温差 逐渐 减 
小 ; 对 比 观察 图 4， 此 时 对 流 换 热 系 数 开始 沿 流动 
方向 逐渐 增加 ， 表 明 换 热 强度 已 开始 恢复 。 
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图 6 起 始 加 热 段 流体 轴 向 速度 与 温度 分 布 
Fig.6 Distribution of axial velocity and temperature in 


the initial heating region. 


进一步 考察 变 物性 及 浮 升力 对 于 起 始 加 热 段 
的 影响 ， 图 7 给 出 保持 流体 密度 、 粘 度 不 变 
恒定 为 入 口 温度 条 件 下 流体 的 物性 )， 以 及 
速度 设 为 零 值 (e=0) 时 计算 内 壁 温 与 实验 数 
对 比 。 如 图 7 所 示 ， 当 流体 密度 保持 恒定 时 
温 未 出 现 局 部 峰值 ， 当 重力 加 速度 设 为 零 值 
内 壁 温 变 化 趋势 与 密度 保持 恒定 时 基本 相 

当 流 体 粘 度 保持 恒定 时 ， 入 口 段 内 壁 温 相 比 
度 恒 定时 有 所 升 高 ， 但 仍 低 于 实际 值 。 由 此 
， 流 体 受 热 后 由 密度 变化 及 重力 引起 的 浮 升 
是 导致 加 热 段 传 热 恶化 的 主要 因素 。 这 一 
与 层 流 混 合 对 流 的 情况 相 类 似 。 在 起 始 加 热 
边界 层 呈 层 流 状态 。 对 于 向 下 流动 ， 浮 升力 
方向 与 流体 流动 方向 相反 ， 导 致 近 壁 区 流 
度 减 小 C21。 其 结果 是 对 流 作用 所 能 传递 的 能 
少 ， 传 热能 力 下 降 ; 壁面 处 速度 梯度 的 减 小 
切 应 力 减 小 ， 推 迟 了 层 流 向 湛 流 的 转 按 ， 二 


pan 


者 综合 作 


]， 最 终 导致 了 传 热 恶 化 的 发 生 。 


此 外 ， 


流体 粘度 的 变化 也 是 造成 入 口 段 传 热 


恶化 的 原因 之 一 。 流 体 受热 后 粘度 减 小 ， 近 壁 


x1 


粘性 层 内 切 应 力 减 小 ， 导 致 了 传 热能 力 的 下 降 。 
但 其 影响 小 于 浮 升力 作用 的 影响 。 
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7. 变 物性 及 浮 升力 对 于 起 始 加 热 段 换 热 的 影响 
Fig.7 Effects of property variations and buoyancy 


forces on heat transfer in the initial heating region. 
+ y 
4 结论 


本 文 使 用 开源 流体 动力 学 程序 库 OpenFOAM 


对 竖 直 圆 管内 超 临 界 压力 RP-3 航空 煤 ; 
热 进行 了 数值 研究 ， 重 点 讨论 了 起 始 加 热 段 的 传 
热 恶 化 现象 。 研 究 主 要 结论 如 下 : 


的 对 流传 


(OD ko SST 模型 无 法 对 起 始 加 热 段 的 传 热 恶 化 


现象 进行 预测 ，k-e-w? 模 型 能 够 定性 描述 起 始 加 
热 段 内 壁 温 的 变化 规律 ， 计 算 结果 与 实验 数据 吻 


合 较 好 。 
(2) 由 流体 密度 变化 及 重力 引起 的 浮 升力 作用 是 


导致 起 始 加 热 段 传 热 恶化 的 重要 原因 。 流 体 粘度 
对 于 起 始 加 热 段 传 热 的 影响 小 于 浮 升力 的 作用 。 
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